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ABSTRACT 



We present a magnetic field mapping of water maser clouds in the star-forming 



region W3 IRS 5, which has been made on the basis of the linear polarization 
VLBI observation. Using the Very Long Baseline Array (VLB A) at 22.2 GHz, 16 
of 61 detected water masers were found to be linearly polarized with polarization 
degrees up to 13 %. Although 10 polarized features were widely distributed in 
the whole W3 IRS 5 water maser region, they had similar position angles of the 
magnetic field vectors (~ 75° from the north). The magnetic field vectors are 
q roughly perpendicular to the spatial alignments of the maser features. They are 

consistent with the hourglass model of the magnetic field, which was previously 
proposed to explain the magnetic field in the whole W3 Main region (r ~ 0.1 pc). 
^ ■ They are, on the other hand, not aligned to the directions of maser feature proper 

motions observed previously. This implies that the W3 IRS 5 magnetic field was 
controlled by a collapse of the W3 Main molecular cloud rather than the outflow 
i ■ originated from W3 IRS 5. 

o ■ 
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1. Introduction 

Birth of massive stars is characterized by a collapse of molecular cloud core when self- 
gravity overwhelms the magnetic pressure in the core as a result of magnetic field dilution. 
In this situation, the core is not supported by the magnetic field pressure in the cloud core 
any more, and called to be magnetically super-critical (e.g., Shu, Adams, & Lizano 1987; 
Feigelson & Montmerle 1999, and references therein). Measurements of the magnetic field 
strengths and directions in star-forming regions have been made in a variety of spatial scales, 
from whole star-forming regions (1-10 pc) down to individual molecular cloud cores (e.g., 
Heiles et al. 1993; McKee et al. 1993; Evans II 1999, and references therein). In the smaller 
scale, the measurements were made mainly with radio interferometers with higher angular 
resolution. However, because of complicated structures of the magnetic fields expected from 
the presence of the dense clumps and outflows, and because of the uncertainty in calibration 
in the polarimetric interferometry, it is more difficult to map the magnetic field accurately 
in this scale. 

Water maser emission is often associated with star-forming regions (SFRs) with energetic 
outflows from young stellar objects (YSOs) (e.g., Elitzur 1992a, b and references therein). 
VLBI observations have revealed that water maser emission consists of clusters of compact 
maser features with a typical size of 1 AU and a velocity width of 1 km s _1 (e.g., Reid & 
Moran 1981). Measurements of Doppler velocities and proper motions of maser features 
have led to the conclusion that the spatial motions of maser features really reflect the true 
three-dimensional gas kinematics (e.g., Gwinn, Moran, & Reid 1992; Torrelles et al. 2001a, 
b). Previous observations have found expansion motions produced by outflows from YSOs 
or expanding HII regions (e.g., Imai et al. 2000, hereafter Paper I, and references therein). 
Shocks, which have velocities exceeding about 20 km s" 1 and are running into high-density 
(n U2 > 3xl0 6 cm" 3 ) magnetized material, successfully explain the water maser emission 
associated with energetic YSO outflows (e.g., Hollenbach & McKee 1989; Elitzur, Hollenbach, 
& McKee 1989, 1992; Hollenbach, Elitzur, & McKee 1993). 

Polarimetric mapping of individual water maser features were made for the first time 
using VLBI data in W51 M by Leppanen, Liljestrom, & Diamond (1998) (hereafter LLD). 
They found that the observed directions of linear polarization were well aligned along a 
stream of water maser features found in W51 M. It should be noted that the observed 
direction of the linear polarization does not straightly indicate the magnetic field direction 
of the cloud because it depends on the angle between the magnetic field and the photon 
propagation direction (toward an observer) (e.g., Goldreich, Keeley, & Kwan 1973; Deguchi 
& Watson 1986b, 1990). In the case of W51 M, the observed directions of linear polarization 
were interpreted to be parallel to the magnetic field vector projected on the sky plane because 



- 3- 



the field vector is expected to be roughly tangential to the line of sight due to coupling 
with the Sagittarius spiral arm. LLD concluded that the alignment of the water maser 
linear polarization was created by shocks caused by the nearby expanding HII region. The 
linear polarization of water masers as a tracer of the interstellar magnetic fields is indirectly 
supported by that the linear polarization directions in a maser region are stable on a timescale 
of several years despite of much shorter time scale of lifetimes and flux variability of maser 
features (Abraham & Vilas Boas 1994; LLD; Horiuchi, Migenes, & Deguchi 2000; Horiuchi 
& Kameya 2000; Horiuchi & Kemaya 2003 in preparation). 

This paper discusses the results of magnetic field mapping of water maser clouds in 
the massive-star forming region W3 IRS 5 revealed by a VLBI observation. W3 IRS 5 is 
considered to be at the very early phase of star formation because of lack of high veloc- 
ity components in the water masers, as suggested by Genzel & Downes (1977). Previous 
observations also supported this interpretation on W3 IRS 5 (e.g., Claussen et al. 1994; 
Tieftrunk et al. 1997; Roberts, Crutcher, & Troland 1997). In a 2"x3" field of W3 IRS 5, 
Claussen et al. (1994) and Tieftrunk et al. (1997) found the distribution of a cluster of 
centimeter-continuum emission sources. Imai et al. (2000) investigated the 3-D kinematics 
of the water masers in W3 IRS 5 and estimated a distance to this region to be 1.83±0.14 kpc. 
In the present paper, the value of 1.8 kpc is adopted as the distance to W3 IRS 5. Imai, 
Deguchi, & Sasao (2002) (hereafter Paper II) investigated the physical condition of the re- 
gion where the water masers are excited. They found that turbulent motions dominate the 
morphology and kinematics on a microscopic scale (down to 0.01 AU). Figure 1 shows the 
distribution and proper motions of the water maser features found in Paper I. 

Troland et al. (1989) and Roberts et al. (1993) (hereafter RCTG) measured the Zeeman 
effect of the HI line in W3 IRS 5. They revealed the magnetic field projected in the line 
of sight on a scale of ~ 0.6 pc, and demonstrated an "hourglass" model of the magnetic 
field. The above scale is roughly equivalent to the core size found in this cloud found in 
the 13 CO line by Hayashi, Kobayashi, & Hasegawa (1989). A pinch of the hourglass was 
found to be located very close to W3 IRS 5 region and to have the magnetic field strength 
higher than 1 mG. Greaves, Murray, & Holland (1994) also made a polarimetric observation 
at sub-millimeter wavelengths and estimated the magnetic field in W3 IRS 5 on a scale of 
0.1 pc; the result supports the hourglass model. Barvainis & Deguchi (1989) found that the 
W3 IRS 5 water masers were highly linearly polarized. Sarma, Troland, & Romney (2001) 
(hereafter STR) and Sarma et al. (2002) (hereafter STCR) measured the Zeeman effect of 
the bright components of the water masers and estimated the magnetic field strength along 
the line of sight to be ~ 30 mG. 

In the present paper, we focus on an issue whether or not the magnetic field measured in 
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water masers traces that found on the larger scale. The cause that determines the shape of 
the large-scale magnetic field is also speculated: either the collapse of the parent molecular 
cloud or the formation of a magneto-hydrodynamical jet. Section 2 describes the polarimetic 
observations using the Very Long Baseline Array (VLBA) and the Very Large Array (VLA) 
of the National Radio Astronomy Observatory (NRAO) 1 and data reduction. The results 
are presented in §3. Section 4 discusses mainly the relation between the measured directions 
of maser linear polarization and possible directions and strength of the magnetic field in the 
observed scale (<5") of W3 IRS 5. 



2. Observations and data reduction 

The polarimetric observation of water masers in W3 IRS 5 was made on 1998 Novem- 
ber 21 using 10 telescopes of the VLBA. The water masers and calibrators (NRAO 150, 
JVAS 0212+735, and 0423-014) were observed for 10 hrs in total. The observed signals 
were recorded in a 4-MHz base band with both right and left circular polarization (RCP, 
LCP) in 2-bit sampling. The recorded data were correlated with the VLBA correlator in 
Socorro, which yielded parallel and cross-hand visibilities with 128 spectral channels. The 
correlated data covered a velocity range of —58 km s _1 < Vlsr< —5 km s" 1 and a velocity 
spacing of 0.42 km s -1 per spectral channel. 

Data calibration and imaging were made using the procedures in AIPS 2 mostly. First, 
the data of the parallel-hand-polarization visibilities were reduced in the usual procedures 
(e.g., Reid 1995; Diamond 1995) and fringe-fitting, self-calibration, and complex bandpass 
solutions were obtained independently in RCP and LCP. Although the water maser emission 
was strong enough to apply the "template method" to calibrate the visibility amplitudes 
(Diamond 1995), the band width of 4 MHz was too narrow to obtain a sufficient number of 
emission-free velocity channels in the template total-power spectrum of the maser emission. 
Instead, we applied tables of the measured system noise temperatures and antenna gains. 
In the self-calibration, we obtained a Stokes / maser map of the velocity component (maser 
spot) at —40.2 km s" 1 and complex gain solutions in every 6 s. 

Second, the instrumental polarization parameters (D-terms) of the telescopes were ob- 
tained using the AIPS RUN script CROSSPOL and the task LPCAL with scans of 0420-014, 
which is well linearly polarized. The D-term solutions were applied to the maser data. After 



1 The NRAO is a facility of the National Science Foundation of USA, operated under a cooperative 
agreement by Associated Universities, Inc. 
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the D-term calibration, the Stokes /, Q, U maps of maser spots were obtained independently. 
The naturally-weighted synthesized beam was 0.68x0.49 mas, with the major axis at a po- 
sition angle (PA) PA = 3°.3. No smoothing in velocity channel was applied. The detection 
limit of about 30 mJy beam -1 at a 5-a noise was obtained in velocity channels without bright 
maser emission. Before making the final image cubes, we looked for the masers by making 
wider maps covering a 2"x5" field with a larger synthesized beam. 

Third, each velocity component (maser spot, see e.g., Gwinn 1994a, b; Paper II) in the 
image cube was fitted by a two-dimensional Gaussian brightness distribution using the AIPS 
task JMFIT. The uncertainty of the estimated position due to thermal noise was 3 microarc- 
seconds (/xas) for a 1-Jy spot. On the other hand, positional accuracy of the brightness peak 
of a maser feature (a cluster of maser spots) depends on the distribution of maser spots in 
the feature and the brightness distributions of the maser spots and limited to be typically 10 
//as. The electric vector position angle (EVPA, see LLD), EVPA = arctan(Q/£/)/2 of each 
maser spot was calculated from the Stokes Q and U parameters. For each maser feature 
(a cluster of maser spots existing together within 1 mas and 1 km s" 1 in position and 
velocity, Paper I and II), the degree of linear polarization and the EVPA were obtained from 
the Stokes parameters (I, Q, V) of the intensity peak in the maser feature. 

The water masers were also observed with the VLA (26 antennas) in C-configuration on 
1998 November 28, one week later the VLB A observation. The water masers and calibrators 
(3C 48, 0420—014, and 3C 138) were scanned for 2 hrs in total. The observed signals were 
recorded in a 3-MHz base band both in right and left circular polarization. The correlated 
data had 64 spectral channels with a velocity spacing of 0.63 km s _1 per spectral channel. 
The D-terms of the antennas were obtained from the scans of 3C 48, which covered a variation 
of parallactic angle over 180°, assuming the intrinsic EVPA of — 72°.5. Then we confirmed 
that the obtained EVPAs of 3C 138 and 0420—024 were consistent with those previously 
obtained by Akujor et al. (1993) and Marscher et al. (2002), respectively. 

The EVPA at VLsr= -42.2 km s' 1 in the maser feature W3 IRS 5:12003 5 (see table 
1) was compared between the VLB A and VLB A data. At this velocity, the maser spot 
was most prominent and its polarimetric parameters were well determined. A correction 
(AEVPA = 39. c 3) was applied to the EVPAs obtained with the VLBA. Taking into account 
the variation of EVPA in the maser feature (< 3°), the ambiguity of the assumed EVPA 
of 3C 48 (~ 5°), and the accuracy of the calibration of the VLA data (uncertainty in R— L 
phase determination, < 0.3°), the accuracy of EVPA was expected to be ~6°. Because the 
VLA observation was made 7 days after the VLBA observation, the estimated EVPAs might 
suffer from the EVPA time variation of the feature W3 IRS 5:12003 5. 
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3. Results 

We detected 61 maser features consisting of 290 maser spots. These numbers are smaller 
by factors of 2-3 and 4-7, respectively, than those in the previous VLBA observations (Paper 
I, see also Paper II) while the sensitivity of the present observation was better than the 
previous ones. The decrease of features was likely due to less activity in maser excitation, 
while the decrease of spots was due to less velocity resolution by a factor of 4 in the present 
observation. 

Table 1 summarizes the parameters of water maser features in W3 IRS 5. Only for 
15 of 61 features, the proper motions were measured previously in Paper I. It was difficult 
to identify the same features between the observations of Paper I and the present work 
because many features were likely to appear and disappear in the crowded regions during 
the epoch separation of 11 month. The highest degrees of polarization was 13% in the 
present observation, which was smaller than those (up to 30%) found in a previous single 
dish observation (Barvainis & Deguchi 1989) and an unpublished VLBI observation with a 
400-km baseline. As shown in figure 2, no correlation was found between a feature intensity 
and a degree of linear polarization. LLD found a systematic change of EVPAs in a single 
maser feature with the full velocity width AV = 1.3 km s _1 . On the other hand, as shown 
in figure 3, EVPAs were almost the same (within 15°) in all of the features with measured 
linear polarization in W3 IRS 5. 

Figure 4 presents the spatial distribution and the linear polarization directions of the 
maser features. The present maser map was easily overlaid on the previous ones in Paper 
I. We found that a cluster of maser features is elongated in the north-west and south-east 
direction and that the radial velocities have been almost the same among previous maser 
maps of Paper I and II and the present one (figure 4c). Very likely this cluster corresponds 
to the cluster of maser features with the measured Zeeman effect (the features in figure 1 of 
STR, see also STCR). As STR mentioned, the cluster was located very close to the radio 
continuum source "c" (Claussen et al. 1994, see in figure 1). It was quite difficult, however, 
to exactly identify the same features from one observing epoch to another because of too 
large time spans between two observations (> 14 months). 

The observed EVPAs were roughly the same among distant maser features; see the 
features in figures 4b and c, within a separation of 3800 AU, as well as adjacent maser features 
within 100 AU. Figure 3 shows that 10 out of 16 polarized features exhibited this tendency. 
The EVPAs were also constant at — 13±5° in the individual features. The exception was 
found in the cluster of features seen in figure 4e, in which EVPAs changed by 90°. Only one 
feature at the south-west corner of the cluster, W3 IRS 5:12003 40, had an EVPA roughly 
equal to those in other clusters. 
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In W51M (LLD) and, EVPAs of masers are well aligned along the direction of a highly 
collimated flow. In W3 IRS 5, on the other hand, EVPAs are not parallel to the direction 
of the whole outflows, which is oriented at PA ~10° (see figure 1 and figure 6 of Paper I). 
In addition, they are not aligned in the direction of the large-scale CO outflow (PA ~ 50° 
, Mitchell, Maillard, & Hasegawa 1991; Hasegawa et al. 1994). Note that the alignment of 
EVPAs observed in the present paper appears more evidently than the alignment of proper 
motions observed in Paper I. This is probably because the random motions are dominant in 
the proper motions. 



4. Discussion 

4.1. Relation between the EVPA and the magnetic field 

The present work has investigated for the first time the magnetic field of the whole 
water maser regions covering a large mapping area (2"x5") and the wide maser velocity 
range (AV = 53 km s _1 ). Despite the proper motions that have quite different PAs among 
the maser features, we found the similar EVPA values (— 13°.7(±3°.7 in median) among 10 
maser features widely distributed in the whole maser region. This position angle should be 
perpendicular or parallel to the magnetic field in the water maser region of W3 IRS 5. 

A maser theory suggests that the directions of the linear polarization of saturated masers 
should be perpendicular (or parallel) to the magnetic field projection on the sky if the 
inclination angle of the magnetic field with respect to the line of sight is \0m\ ^ 55° (or 
\0m\ ^ 55°) (Goldreich, Keeley, & Kwan 1973; Deguchi & Watson 1986b). This inclination 
angle was roughly estimated from the observable parameter Q/I (see figure 2 of Deguchi & 
Watson 1986a). We estimated the inclination angle to be 45° < M < 60° for the observed 
value \Q/I\ < 0.15, which is close to the critical inclination angle mentioned above (~ 55° 
). We expect that the feature at the south-west corner of the cluster, W3 IRS 5:12003 40, 
exhibits a "flip" of the linear polarization direction by ~ 90° because the magnetic field 
has an inclination angle close to the critical value. However, in other features, we favor 
\0m\ ^ 55° because the magnetic field on the larger scale is almost tangential to the sky 
plane as mentioned above. Thus, the magnetic field projection was estimated to have a 
position angle PA ~ 76°. In contrast, the directions of linear polarization in the W51 M 
water maser are considered to be parallel to the magnetic field (LLD). 

Note that the linear polarization is affected more or less by Faraday rotation due to free 
electron in the interstellar medium along the line of sight. We evaluated Faraday rotation 
for the maser linear polarization at 1.35-cm wave length, $, using the following qeuation, 
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$ ~ 1.48 x 10 2 Dn c H\\, (1) 

where D (pc) is the size of an interstellar cloud, n c (cm -3 ) the number density of free 
electron, and H\\ (G) the magnetic field strength along the line of sight. For the cases of a 
maser cloud (D ~ 1CT 4 pc, n c ~ 1 cm -3 , H\\ ~ 30 mG), a molecular cloud (D ~ 0.1 pc, 
n c ~ 1 cm" 3 , if || ~ 100 //G), and the interstellar space between the W3 region and the Sun 
(D ~ 1.8 kpc, n c ~ 10~ 2 cm" 3 , H\\ ~ 10 fiG), we obtain $ ~ 4.4xl0~ 4 rad, 1.5xl0- 3 rad, 
and 2.7xl0~ 2 rad, respectively, which are negligible effects. Although the hyper-compact HII 
regions found in the W3 IRS 5 region (D ~ 2xl0~ 3 pc, n c ~ 1.3 xlO 5 cm -3 , H\\ ~ 30 mG, 
Claussen et al. 1994) might introduce the strongest Faraday rotation ($ ~ 120 rad), they 
are not located in front or behind the water maser features with linear polarization. 



4.2. The magnetic field on the arcsecond-scale in W3 IRS 5 

Troland et al. (1989), RCTG, and Greaves, Murray, & Holland (1994) found a gradient 
of the magnetic field along the line of sight and proposed "hourglass" models for the W3 IRS 5 
magnetic field on a large scale (1'). In these models, an inclination angle of the hourglass 
major axis 9 ~ 2° was adopted with respect to the sky plane on the basis of the angular 
distribution of the magnetic field strength along the line of sight, which was obtained by 
RCTG. A position angle of the major axis, PA ~ 50°, and a FWHM width of the pinch, 
FWHM = 10", were adopted on the basis of the spatial distribution of the EVPAs, which 
was obtained by Greaves, Murray, & Holland (1994) with an angular resolution of 13". 5. 

In this paper, we propose an hourglass model of the magnetic field shown in figure 5. 
This model is made to explain the magnetic field on the small scale (< 10"), and consistent 
the large scale magnetic field (~ 1') observed by RCTG and Greaves, Murray, & Holland 
(1994). The modeled field is pinched- in at the peak of the 800-/zmcontinuum emission. 
Ladd et al. (1993) observed the 450-/zm/800-/zmcontinuum emission with a 7"-beam (a 
pointing accuracy of 2"), and found that the peak of the emission has a position offset 
of ~ —2". 5 in R.A. and ~ +0".5 in decl. with respect to the center of the maser cluster. 
Although the continuum emission peak should coincide with the 13 CO and C 18 emission 
peaks (Oldham et al. 1994; Tieftrunk et al. 1995; Roberts, Crutcher, & Troland 1997), the 
observed position is somewhat uncertain because of the large beams. It is possible that the 
C 13 emission peak observed by Roberts, Crutcher, & Troland (1997) traces a segment of the 
outflow while the continuum emission peak traces the densest part of the parent molecular 
cloud. Tieftrunk et al. (1995) found four gas clumps traced by C 34 S J =3-2 and J =5-4 
emission, but none of the emission peaks coincides with the pinch center. 
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The magnetic field projection on the sky (PA ~ 76°) mentioned in the previous sub- 
section is well consistent with that in the proposed magnetic field model (see figure 5b). 
Taking into account that the magnetic field vector is oriented toward the observer (a negative 
value of the measured Zeeman effect, STR; STCR), the maser region is likely located in front 
of the major axis of the hourglass (see figure 5c, d). The inclination of the magnetic field 
to the line of sight, however, cannot be completely consistent with the model, in which the 
inclination significantly varies from a point to another in the whole maser region. 

The alignment of the magnetic field in the coverage of about l"x4" in the water maser 
region provides a constraint on the pinch ratio of the hourglass model around A/B ~ 0.14, 
as well as on the pinch width around FWHM ~ 10") (definitions of A and B are given in 
the caption of figure 5). As discussed by RCTG, one obtains a magnetic field strength of 
~1.3 mG at the pinch center when assuming a strength of the total undisturbed magnetic 
field of 25 /iG and considering the r~ 2 dependence of the magnetic field. This value is one 
order of magnitude smaller than that obtained by STR and STCR (~ 30 mG) but one order 
of magnitude larger than that estimated by RCTG. Crutcher (1999) found an empirical 
relation of a magnetic field strength, B oc n 0A7 , over a large range of gas densities n up to 
10 7 cm" 3 . Therefore, it is expected that the field strength found in the water maser cloud 
(n > 10 8 -10 9 cm" 3 ) is consistent with the B-n law (c.f. STR). STCR also pointed out 
that the magnetic field in the interstellar space (preshock region) should be enhanced in the 
maser region (postshock region) by a factor of ~ 20. 

4.3. The magnetodynamically super-critical collapse of the parent molecular 

cloud of W3 IRS 5 

The origin of the magnetic field pinch is still obscured. The direction of the whole 
magnetic field is almost parallel to a molecular outflow found in W3 A (Mitchell, Maillard, 
& Hasegawa 1991; Hasegawa et al. 1994). Momose et al. (2001) pointed out that a magnetic 
field would control the direction of a molecular outflow rather than that the magnetic field 
is influenced by the outflow. In W3 IRS 5, however, although the magnetic field seems to 
along the molecular outflow on a larger scale (~ 1'), it may not necessarily be parallel to 
the outflows on a smaller scale (~ 5") as seen in the water masers. We evaluated the energy 
densities of the magnetic field and of the outflow of W3 IRS 5 using the same equations 
described in Momose et al. (2001). The critical magnetic field strength, f? critical , at which 
the field has the same energy density as that of the cylindrical outflow, is expressed by 
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where E aovi , r and I are the kinetic energy, the radius, and the total length of the cylindrical 
outflow, respectively. Adopting the total outflow energy estimated by Hasegawa et al. (1994), 
E ao „ ~ 1.2 xlO 46 erg, and the radius and the length of the outflow estimated by Mitchell, 
Maillard, & Hasegawa (1991), r ~ 4.1xl0 17 cm and I ~ 2.5xl0 18 cm, respectively (assuming 
the distance to W3 IRS 5 of 1.8 kpc here), we obtained -B crit icai — 0.48 mG, which is larger 
than the maximum strength observed in the HI Zeeman effect by RCTG (< 0.22 mG). To 
obtain -B critica i for the water maser cloud, we converted equation (2) to, 



where p mascr > 3.3xl0~ 14 g cm -3 (corresponding to n H2 = 10 8 cm -3 ) is the gas density of 
the maser region necessary for maser excitation, t> maser is the outflow velocity. Adopting the 
outflow velocity f maser ~ 15 km s -1 , we obtained -B critica i ~ 250 mG. Thus, the magnetic field 
strength estimated on the basis of the Zeeman splitting (~ 30 mG, STR) is one order of 
magnitude lower than the critical strength. This suggests that the magnetic cannot control 
the outflow dynamics. 

This evaluation is consistent with that made for the W51M region by LLD. On the 
other hand, there is no observational evidence that the outflow controls the magnetic field 
structure. Momose et al. (2001) suggested that the magnetic field is aligned more in the 
younger phase of star-formation than in the more evolved one so that the disturbance of the 
former magnetic field due to the small-scale turbulence is not significant yet. This suggestion 
supports that W3 IRS 5 should be at the very early stage of star formation. Oldham et al. 
(1994) and Tieftrunk et al. (1995) estimated the mass of the molecular cloud associated with 
W3 IRS 5 to be 800-900 M Q . This value is still larger than the magnetodynamically-critical 
mass of 510 M Q for a cloud with a radius of 9000 AU (corresponding to 5" at 1.8 kpc) and 
a magnetic field strength of 30 mG. This indicates that the magnetic field of the W3 IRS 5 
cannot support the gravity of the molecular cloud, in other words, that a part of the W3 IRS 5 
cloud is collapsing under the magnetodynamically super-critical condition (Shu, Adams, & 
Lizano 1987). 

Note that the above evaluation is still consistent with the suggestion that the magnetic 
field should control the water maser excitation (Elitzur, Hollenbach, & McKee 1989, 1992; 
STCR). To consider the maser excitation, the velocity dispersion of the postshock region 
where the masers are excited, a v < 0.6 km s -1 , is inserted instead of f maser in equation 3 
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(Paper I, II; STCR). Then one can estimate -B masor ~ 10 mG, smaller than or equal to the 
observed magnetic field strength. 

It is still unclear what kind of object exists in the pinch center of the magnetic field. 
Tieftrunk et al. (1995) found four molecular gas clumps in W3 IRS 5, but none of the clumps 
is located at the pinch center. It is likely that these clumps were formed in the common 
parent molecular cloud and are responsible for the observed magnetic field. The magnetic 
field which controlled the multiple clumps supports the mild pinch ratio A/B ~ 0.14 rather 
than the extremely small ratio A/B ~ 0.04 speculated by RCTG. 

5. Summary 

Polarimetric VLBI observations of water masers in W3 IRS 5 have revealed that observed 
EVPAs of water masers features are well aligned in the whole maser region except for a small 
number of maser features. The directions of the magnetic field projected on the sky in the 
maser region are estimated to be perpendicular to the directions of the linear polarization 
of the masers {PA ~ 76°). The magnetic field estimated supports the hourglass magnetic 
field model proposed previously. Thus, water maser linear polarization may reveal in detail 
the magnetic fields of the interstellar medium in star-forming regions. The W3 IRS 5 water 
maser region has a position offset from the pinch center of the hourglass and its magnetic 
field is unlikely to control the outflow dynamics. All of these observational facts suggest 
that the W3 IRS 5 cloud is collapsing under magnetically super-critical condition and that 
the outflows are not yet developed enough to completely disturb the magnetic field in this 
cloud. The hourglass structure of the magnetic field is suggested to be formed as a result of 
the collapse of the parent molecular cloud that is magnetically super-critical. 
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Table 1: 


Parameters of the water 


maser 


features toward W3 IRS 5 






Mascr feature 1 


Offset 2 




LSR velocity (km s 1 ) 


Peak Intensity 


pi 


EVPA 


Mascr feature' 1 


(W3 IRS 5 


12003) 


R.A. (mas) decl. (mas) 


Peak 


Width 3 


(Jy beam - 1 ) 


(%) 


(dog) 


(W3 IRS 5: 12000) 


1 




-1.50 


3.33 


-56.0 


1.7 


0.04 








2 




1075.72 


1016.29 


-43.1 


3.8 


18.41 


11.0 


-12.1 




3 




-389.93 


-661.34 


-42.9 


1.7 


7.78 








4 




0.08 


-0.99 


-42.6 


1.3 


6.95 








5 




-1.49 


3.36 


-42.1 


3.8 


111.05 


0.8 


-15.1 




6 




-388.32 


-665.56 


-42.0 


2.1 


9.60 


0.7 


19.8 




7 

8 




-388.08 
-431.85 


-671.67 
-620.46 


-41.7 
-41.7 


0.8 
0.8 


3.63 
1.48 








9 




-387.87 


-667.89 


-41.3 


0.4 


0.50 








10 




-0.01 


-0.09 


-41.3 


1.3 


26.10 








11 




-18.42 


0.16 


-40.7 


1.7 


8.95 








12 




-17.26 


-1.79 


-40.5 


1.3 


2.63 








13 




-1.37 


4.38 


-40.5 


1.7 


8.69 








14 




0.00 


0.00 


-40.5 


5.1 


79.89 


0.8 


-14.9 




15 




-1.66 


5.47 


-40.3 


1.3 


8.91 








16 




-0.37 


-11.64 


-40.0 


2.5 


20.81 








17 




-3.92 


5.23 


-39.6 


0.4 


5.05 








18 




-387.50 


-670.38 


-39.3 


4.2 


9.99 


0.8 


-17.8 




19 




-386.67 


-677.37 


-38.9 


3.0 


7.25 








20 




-4.28 


8.23 


-38.8 


1.3 


5.63 








21 




-3.78 


7.50 


-38.8 


0.8 


3.91 








22 




-6.27 


8.67 


-38.2 


3.8 


26.49 


0.7 


-11.8 




23 




9.72 


-16.39 


-37.8 


2.1 


11.11 


12.8 


-11.4 




24 




-386.73 


-676.74 


-37.7 


1.3 


0.29 








25 




-7.97 


13.24 


-37.5 


1.7 


1.88 








26 




-9.20 


14.11 


-36.9 


1.7 


2.14 








27 




786.55 


1792.69 


-36.7 


2.1 


10.01 


1.1 


-15.5 


64 


28 




797.51 


1750.03 


-36.2 


0.8 


0.56 






65 


29 




783.13 


1802.02 


-35.8 


3.8 


14.81 


0.9 


-18.0 


66 


30 




787.03 


1794.17 


-35.5 


5.1 


91.10 


0.4 


-12.5 


63 


31 




328.54 


-2001.74 


-35.0 


1.7 


2.24 






71 


32 




796.47 


1750.38 


-35.0 


3.0 


0.53 








33 




-416.12 


-600.11 


-34.6 


4.2 


11.55 








34 




-1.46 


3.30 


-34.6 


0.4 


0.11 








35 




-0.02 


-0.02 


-34.6 


0.4 


0.06 








36 




-6.21 


8.51 


-34.6 


0.4 


0.08 








37 




794.01 


1751.23 


-34.5 


0.4 


2.10 








38 




794.65 


1750.92 


-34.1 


0.4 


0.23 






79 


39 




295.92 


-2025.00 


-33.9 


1.7 


0.94 






73 


40 




294.72 


-2025.82 


-33.8 


2.5 


4.43 


5.4 


-5.0 


76 


41 




293.69 


-2027.11 


-33.7 


0.8 


1.06 








42 




331.95 


-2004.59 


-33.7 


1.3 


6.95 




89.6 




43 




325.76 


-1999.60 


-33.6 


4.6 


21.89 


7.3 


30.4 


82 


44 




332.52 


-2004.95 


-33.3 


0.8 


4.25 


12.9 


-79.2 




45 




-415.21 


-602.42 


-33.2 


3.0 


1.00 








46 




334.25 


-2006.04 


-33.1 


3.0 


10.90 


8.4 


-59.1 


83 


47 




-415.45 


-604.25 


-33.0 


3.4 


3.84 








48 




324.56 


-1998.90 


-32.9 


5.9 


13.12 








49 




324.99 


-1998.91 


-32.4 


2.1 


6.56 


1.5 


-62.5 




50 




795.87 


1750.48 


-32.3 


2.1 


1.93 








51 




326.65 


-1999.60 


-32.0 


1.3 


0.19 






77 


52 




331.83 


-2004.52 


-31.2 


0.4 


0.05 








53 




249.87 


-2079.93 


-30.4 


4.6 


4.02 






91 


54 




228.02 


-2083.62 


-29.9 


1.7 


1.05 






94 


55 




230.35 


-2083.38 


-29.9 


0.8 


0.44 






98 


56 




229.26 


-2083.48 


-29.9 


0.4 


0.01 








57 




227.05 


-2083.52 


-29.9 


1.7 


1.09 








58 




325.72 


-1999.57 


-29.6 


2.1 


0.06 








59 

60 




795.39 
-503.83 


1750.80 
-2377.04 


-29.5 
-23.8 


1.3 
3.0 


0.13 
0.58 








61 




-205.57 


-1830.29 


-6.8 


2.1 


1.32 









1 Water maser features detected toward W3 IRS 5. The feature is designated as W3 IRS 5:12003 N, where N is the ordinal source number given 
in this column (12003 stands for sources found by Imai ct al. and listed in 2003). 

2 R.clative value with respect to the location of the position-reference maser feature: W3 IR.S 5:12003 14- 
^ A channel velocity spacing multiplied by a number of detected maser spots in the feature. 

4 Fraction of linear polarization intensity to total intensity. 

5 The feature identical to that found in Paper I, W3 IRS 5:12000 M, where M is the ordinal source number given in this column. 
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Relative R.A. offset (mas) 

Fig. 1. — Water maser features and compact radio continuum sources in W3 IRS 5. Arrows 
indicate locations and proper motions of maser features, which were revealed by Paper I. The 
arrow length indicates the magnitude of the proper motion per year (x325), respectively. 
The radio continuum sources, which were mapped by Claussen et al. (1994), are shown as 
contours. The designations the continuum sources are same as those in Tieftrunk et al. 
(1997). Two plus marks A and B indicate the locations of driving sources of the outflows in 
W3 IRS 5, which were estimated by Paper I. 



14 



£ 12- 
g 10 



8 - 

O 



"s. 6 
S 

s 



o 
a. 

s 

— 



4 



2 - 



20 40 60 80 100 120 

I , (Jybeam 1 ) 

peak ^ ' 



Fig. 2. — Distribution of polarization degrees of maser features against their peak intensities. 
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Fig. 3. — Distribution of EVPAs of maser spots against their Doppler velocities. A filled 
circle shows the maser spot detected its linear polarization. Maser spots in the same maser 
feature are connected each other with a thick line. 
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Fig. 4. — Map of linearly-polarized maser features in W3 IRS 5. (a): Locations of the 
maser features (filled circle) detected on 1998 November 21. Two plus marks are the same 
as those in figure 1. The locations of the maser features were estimated with respect to 
the driving sources of the outflows by assuming the stable distribution of clusters of maser 
features among observations made in the previous (Paper I) and the present works. The 
location uncertainty of the A outflow will be within 100 mas. (b)-(e): Details of the feature 
distribution and the directions of linear polarization of the features. A length and a position 
angle of a bar indicate a logarithmic scale of the degree of linear polarization and an EVPA, 
respectively. 
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Fig. 5. — Hourglass model of the magnetic field in the W3 IRS 5 region. The angle 9 ~ 2° is 
the angle of inclination with respect to the sky plane, ~ 50° the position angle, A/B the 
pinch ratio (A = 5", B = 35") and FWHM = 10" the FWHM width of the the pinch at 
the location of the water masers. The all of the definitions are the same as those described 
in figure 7 of Robert et al. (1993). A square (10"xl0") in the overview panel indicates the 
area zoomed up and described in other three sub-panels. Thick parts of the magnetic field 
lines in the zoomed-up front view sub-panel indicate the locations where the field line have 
position angles 65° < PA M < 85°. The location and the area of the W3 IRS 5 water maser 
region are indicated by a gray ellipse. The area center is located at 2". 5, — 0".5, and — 2".0 
in R.A., decl., and line-of-sight offsets, respectively. The orientation of the ellipse major axis 
indicates the direction of the outflow found in the maser region (see figure 6 and 7 of Paper 
I). Thick parts of the magnetic field lines in the zoomed-up top and east-side views indicate 
the locations where the field line have inclinations 45° < \9 M \ < 60°. 



